
use data, clear
set more off
set scheme s1color


gen expelleed=expellee
gen expelleea=expellee

local controls i.sex i.state i.F4 i.F13 i.age10 


estimates clear
reg donate_st expelleed `controls' if t2==1&above==0, ro
eststo m_pess_don
reg logamount expelleea `controls' if t2==1&above==0, ro
eststo m_pess_am
reg pc_att expellee  `controls' if t2==1&above==0, ro
eststo m_pess_att

reg donate_st expelleed   `controls' if t2==1&above==1, ro
eststo m_opt_don
reg logamount expelleea   `controls' if t2==1&above==1, ro
eststo m_opt_am
reg pc_att expellee  `controls' if t2==1&above==1, ro
eststo m_opt_att


coefplot (m_pess_don, aseq("Donate") keep(expelleed) mcolor(blue) ciopts(lcolor(blue blue)) label("Pessimistic")) ///
(m_opt_don, offset(-0.15) aseq("Donate") keep(expelleed) mcolor(red) ciopts(lcolor(red red)) label("Optimistic")) ///
(m_pess_am, aseq("Log Amount") keep(expelleea) mcolor(blue) ciopts(lcolor(blue blue)) nokey) ///
(m_opt_am, offset(-0.15) aseq("Log Amount") keep(expelleea) mcolor(red) ciopts(lcolor(red red)) nokey) ///
(m_pess_att, aseq("Attitudes") mcolor(blue) keep(expellee) ciopts(lcolor(blue blue)) nokey) ///
(m_opt_att, aseq("Attitudes") offset(-0.15)  mcolor(red) keep(expellee) ciopts(lcolor(red red)) nokey), ci(95) xline(0, lcolor(gs10)) ///
ysca(titlegap(2)) msize(small) ///
xtitle("Difference expellee - non-expellee", size(small)) ytitle("", size(small)) ///
xlabel(, nogrid labsize(small)) ylabel(, nogrid labsize(small)) xsca(titlegap(2)) ///
grid(none) swapnames legend(size(small) cols(1)) title("Without information", size(small)) name(noinfo)



estimates clear
reg donate_st expelleed `controls' if t3==1&above==0
eststo m_pess_don
reg logamount expelleea `controls' if t3==1&above==0
eststo m_pess_am
reg pc_att expellee  `controls' if t3==1&above==0
eststo m_pess_att

reg donate_st expelleed   `controls' if t3==1&above==1
eststo m_opt_don
reg logamount expelleea   `controls' if t3==1&above==1
eststo m_opt_am
reg pc_att expellee  `controls' if t3==1&above==1
eststo m_opt_att


coefplot (m_pess_don, aseq("Donate") keep(expelleed) mcolor(blue) ciopts(lcolor(blue blue)) label("Pessimistic")) ///
(m_opt_don, offset(-0.15) aseq("Donate") keep(expelleed) mcolor(red) ciopts(lcolor(red red)) label("Optimistic")) ///
(m_pess_am, aseq("Log Amount") keep(expelleea) mcolor(blue) ciopts(lcolor(blue blue)) nokey) ///
(m_opt_am, offset(-0.15) aseq("Log Amount") keep(expelleea) mcolor(red) ciopts(lcolor(red red)) nokey) ///
(m_pess_att, aseq("Attitudes") mcolor(blue) keep(expellee) ciopts(lcolor(blue blue)) nokey) ///
(m_opt_att, aseq("Attitudes") offset(-0.15)  mcolor(red) keep(expellee) ciopts(lcolor(red red)) nokey), ci(95) xline(0, lcolor(gs10)) ///
ysca(titlegap(2)) msize(small) ///
xtitle("Difference expellee - non-expellee", size(small)) ytitle("", size(small)) ///
xlabel(, nogrid labsize(small)) ylabel(, nogrid labsize(small)) xsca(titlegap(2)) ///
grid(none) swapnames legend(size(small) cols(1)) title("With information", size(small)) name(info)


graph combine noinfo info, xcommon iscale(1)

